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Abstract. Semiclassical techniques have proven to be a very powerful method to extract 
physical effects from different quantum theories. Therefore, it is expected that in the near 
future they will play a very prominent role in the context of quantum gravity. In this work we 
develop systematic tools to derive semiclassical approximations for any quantum theory with 
one degree of freedom. In our approach, the wave function is decomposed in terms of an infinite 
set of moments, which encode the complete quantum information of the system. Semiclassical 
regimes can then be properly described by truncation of this infinite system. The use of efficient 
computer algebra tools allows us to compute the equations of motion up to a very high order. 
In this way, we can study very precisely the quantum back reaction of the system as well as 
the convergence of the method with the considered order. Finally, these tools are applied to 
the particular case of a homogeneous universe filled with a massless scalar field and positive 
cosmological constant, which provide interesting physical results. 



1. Introduction 

Nowadays, one of the main goals of the research in the context of quantum gravity is to extract 
physical predictions from the different theories proposed in the literature. This is a crucial step 
since, in the best scenario, we will be able to make contact with the experiments in order to 
verify or falsify these theories. As in any other quantum theory, semiclassical approximations 
are expected to play a very prominent role in this respect. On the one hand, they will allow us 
to extract new physical phenomena in their validity regime, where quantum effects are relevant 
but not dominant. On the other hand, even if deep quantum phases can be conceptually very 
interesting, their direct experimental observation is not expected in a near future due to technical 
difficulties. Thus, in order to obtain potentially observable effects, one should rather focus on 
the semiclassical region. In the particular case of quantum cosmology, these approximations 
should provide consistency with the large-scale behavior of the universe and, probably, some 
potential low energy observation, for instance through the cosmic microwave background. 

In this regard, the main available tool in quantum gravity is that of effective equations pT*, "2]. 
This approach describes the evolution of the expectation values of the quantum observables of 
the system as well as its fluctuations, correlations and higher-order moments. In the general case 
they form an infinite system of couple differential equations. However in semiclassical regimes 



it can be truncated and provides, to a very good approximation, the evolution of the quantum 
system. 

Therefore, the main goal of this work is to derive systematic semiclassical approximations for 
any quantum system with one degree of freedom. In order to do so, we will first analyze the 
general formalism and then construct efficient computer algebra tools. Finally, the developed 
machinery will be applied to the particular case of a homogeneous universe with a matter content 
of a massless scalar field and positive cosmological constant. 

2. General formalism 

Let us assume we have a classical system with one degree of freedom which is thus completely 
described by a couple of conjugate variables (y,P). Once the quantization of this system is 
performed promoting those variables to operators, the observable information must be extracted 
through expectation values. Hence, we define the infinite set of moments, 

G»''':=((P-Pr(F-F)%,yi, (1) 

with V := {V) and P := {P), which will encode the complete quantum information of the 
system. These moments are purely quantum variables that would vanish in a classical regime. 

The evolution of these variables is then given by the effective Hamiltonian, defined as the 
expectation value of the quantum Hamiltonian H. By Taylor expansion, it can be written as 

H^,{V, P, G"^'^) := (i7)we,i = H{V, P) + E E (2) 

where H is the classical Hamiltonian. Making use of the definition {{X), {Y)} = —ih{[X,Y]), 
for any two operators X and Y, the equations of motion generated by this effective Hamiltonian 
are guaranteed to be equivalent to the Schrodinger flow of states or the Heisenberg flow of 
operators. 

As has been made explicit, the effective Hamiltonian is a state-dependent object and hence 
a function of expectation values and all the moments G"'^. Therefore, the computation of the 
equations of motion for different variables implies calculating the Poisson bracket between two 
generic moments. This bracket is given by, 

n ^ 2 

where the sum over n runs over all odd numbers from 1 to Min(a + c, 6 + d, o + 6, c + d) and the 
coefficients i^^ftcd depend on all its indices. This formula was already presented in the literature 
[T], but a typo remain discovered. Due to the importance of this formula to systematize the 
computation of the effective equations, in Ref. we have corrected it and provided a complete 
proof. 

3. Algebraic implementation 

In order to automatize the computation of Poisson brackets between two generic moments 
{G"'*, C^'*^}, we have developed two different Mathematica codes. The first code works iteratively. 
That is, it begins with the definition of the moments given by ([1]) and then uses the generic 
relation {{X),{Y)} = —ih{[X,Y]), as well as properties of the Poisson brackets and of the 
commutator, until all terms are reduced to simple commutator between isolated operators V 
and P. This is a long process and, in fact, the time the computer requires to compute a given 
Poisson bracket {G'^'*, G'^''^} by this method increases exponentially with the considered order 
{a + b + c + d). 



On the other hand, the second code makes use of the formula ([3]) and it is much more 
efficient. For instance, the timing to compute a 40th-order bracket hke {G^^'^, G^'^^} is less 
than a thousandth of a second, which clearly shows the practical importance of equation ([3]). 
This generic code is one of the main results of our research and it allows us to obtain the 
equations of motion to a very high-order for any quantum system with one degree of freedom. 
In the following section we will develop one example in the context of quantum cosmology. 



4. An example: massless scalar field with positive cosmological constant 

The Friedmann equation for a homogeneous universe with a massless scalar field cj) and a positive 
cosmological constant A reads, 

where a is the scale factor, p^f, the conjugate momentum of the scalar field and derivatives are 
respect to cosmological time. We define the canonical gravitational variables V := 2a^/3 and 
P := a' /a, and choose the scalar field as the internal time. In this way, at classical level, the 
momentum turns out to be the physical Hamiltonian, 



H:=P4, = -VVP^-A. (5) 

In order to obtain the effective Hamiltonian that will describe the dynamics of all moments, 
the Taylor expansion ^ is applied 

oo 

:= i7 + ^[ya„(P,A)G"'° + /3„(P,A)G"-i'i]. (6) 

n=2 

Here a„ and Pn are coefficients that only depend on the order n, the cosmological constant A 
and the momentum P. Note that, due to the linear dependence of the classical Hamiltonian on 
the variable V, only moments of the form G"'" and G°'^ appear on the definition of the effective 
Hamiltonian. 

In order to extract physical information from the system, we have obtained the equations of 
motion and solve them numerically at every order up to ten. As initial state we have chosen an 
unsqueezed Gaussian pulse in the volume V. The corresponding moments G"''' are non-vanishing 
only if both a and b are even numbers. The initial values of the expectation values V and P, 
as well as the width of the Gaussian are chosen so that initially the state is picked, that is, the 
relative fluctuations AV/V and AP/P are small. Finally, three different values have been taken 
for the cosmological constant: small, intermediate and large. 

The behavior of this system with vanishing cosmological constant (A = 0) is particularly 
simple. In this case, the effective Hamiltonian is only composed by a flnite amount of terms and 
all equations of motion decouple. Therefore there is no back reaction and the evolution of all 
variables is simply exponential V ^ e'^, P ~ e~'^, and G""'^ ~ eib-a)(l)^ This harmonic system is 
used as a reference to check that the quantum back reaction becomes stronger when considering 
larger values for A. On the other hand, the negative cosmological case has been analyzed with 
these effective techniques in [3] and with a complete numerical solution of the wave equation in 
[5]. However, for A < 0, the volume is bounded by recollapse and the moments do not grow 
for large volumes. These two features are not present for A > case, making this analysis of 
high-order moments necessary. 

Regarding the evolution of the expectation values V and P, we would like to point out two 
main features we have observed. Classically, the volume V{(j)) of the universe is an increasing 
function and reaches infinity at a finite value of time (j) = '?^div Our analysis shows that the 



quantum back reaction enhances this divergence. That is, at higher order the divergence of 
the volume is reached faster. The second feature is related to the convergence. When we 
truncate the system at a given order, we implicitly assume that the higher-order moments can 
not change significantly the observed behavior. In order to test this convergence, we have 
defined AVn{(p) := 1 — Vn{(p)/Vn+i{(j)), where Vn{4>) is the trajectory of the expectation value of 
the volume at order n. This object is time dependent, so its behavior will change with (/>. We 
have studied three different values of (j) and in all the considered cases, even for very late times, 
and for all values of A, we have obtained an exponential convergence of the system with the 
order n. This is a very important result and shows the reliability of our method even near the 
divergence of the volume. 

As of the behavior of the moments, we observe that they depart from their harmonic behavior 
(A = case) as the value of the cosmological constant is increased. As expected, several moments 
grow and become dominant, providing thus strong quantum back reaction. In all cases, even 
for small values of A, as soon as the evolution is started, the initially vanishing moments are 
immediately excited. In this way, the initial Gaussian shape is broken. Remarkably, after a 
rapid initial increase of these moments, the state adapts itself to the evolution and moments 
change less severely. One of the open questions in this respect is whether this natural state is 
independent of the initial conditions. We can not give any definitive answer without performing 
further numerical analysis but, due to the rapid adaptation to the evolution, the consequences 
of assuming an initial Gaussian state might, after all, be not so special. 

5. Conclusions 

We have presented a systematic method to analyze the effective behavior of any quantum system 
with one degree of freedom. The key ingredients of this method are the generic formula for the 
Poisson bracket between any two moments ^ and an extensive use of computer algebra tools. 
As an example, these tools have been applied to the particular case of a homogeneous universe 
with a matter content of a massless scalar field and positive cosmological constant. The main 
results we have obtained is that the quantum back reaction enhances the classical divergence of 
the volume. In addition, we have observed that the initial Gaussian state is immediately adapted 
to another state better suited to the evolution. And remarkably, exponential convergence of the 
results with the considered order has been obtained even for late times (near the divergence of 
the volume), which validates the applied methods. 
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